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In this summary we give a concise overview of the experimental and theoretical re- 
sults, which were presented during the QCD and Hadronic Final State Working Group 
sessions at the DIS 2007 workshop. 

1 Introduction 

Final states provide a unique source of information to study elementary particle interactions. 
They make it possible to test our general understanding of QCD, explore QCD evolution 
in limiting regimes and also provide a solid reference for new physics searches. It is beyond 
question that precise input from the Standard Model (SM) is needed for many new physics 
searches that could be carried out at the upcoming LHC, most notably precise values of a s 
and of the parton density functions (PDF) are needed. The role of HERA and the Tevatron 
in this respect is indisputable. In this document we summarize the current experimental 
and theoretical efforts in this direction presented at our Working Group [T] . 

2 Theoretical progress 

2.1 Higher orders and formal developments 

One of the highlights of our Working Group was the presentation of the full NNLO calcu- 
lation of e + e _ — > 3 jets by Aude Gchrmann [2]- The main motivation for this calculation 
is that the current error on a s (Mz) = 0.121 ± 0.001(exp) ± 0.005(th) is dominated by 
theory and that this uncertainty largely comes from missing higher orders. A reduction 
of this error can be achieved by computing NNLO 3-jet event shapes and jet rates. This 
means computing 2-loop matrix elements with 3 partons in the final state, one-loop matrix 
elements with 4 partons and tree-level five parton final states. Additionally one needs to 
develop a formalism to cancel divergences before clustering partons into jets, this is done 
with a NNLO antenna subtraction [3] . After implementing results in an extended version of 
the program EERAD2, one can obtain a thrust distribution at NNLO for all possible color 
structures. The results presented show that large cancellations between the various color 
configurations occur, high precision results are therefore needed. This work is currently in 
progress, together with an implementation of other event shapes and 3-jet rates. This will 
allow a first NNLO determination of a s from event shape data. 

An essential ingredient in performing fixed order calculations within any quantum field 
theory is a method to regulate the ultraviolet (UV) divergences appearing in perturbative 
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calculations. One of the most widely used schemes is dimensional regularization (DREG), 
where one extends the number of dimensions to D = 4 — 2e. It has long been known however 
that dimensional regularization breaks supersymmetry (SUSY), one of the promising new 
physics scenarios which might be discovered at the LHC, because the number of degrees 
of freedom for fermions and bosons is different. Precision calculations within SUSY are 
most conveniently done in a scheme which preserves the boson-fermion symmetry. An 
alternative to DREG is provided by dimensional reduction (DRED), where one splits the 
four-dimensional gluon field in a Z?-dimensional component and a so called e-scalar, a scalar 
whose coupling to fermions vanishes in four dimensions (evanescent coupling). 

One attractive feature of SUSY is the unification of couplings close to the grand uni- 
fication (GUT) scale. If one wants to test unification within SUSY one needs to relate a 

five-flavour MS coupling at the electroweak scale, c4 MS (Mz), where all heavy degrees have 

been integrated out (decoupled), to a full DR coupling at the GUT scale, ai full ^ DR (MGUT)- 
To obtain a next-to-next-to-leading logarithmic (NNLL) result one needs three-loop QCD 
running of the MS coupling from Mz up to Mgugy, the typical mass scale of SUSY parti- 
cles, two-loop decoupling and three-loop running in SUSY from MgusY up to Mqut- Robert 
Harlander reported on this involved calculation [4]. The result turns out to be almost in- 
dependent of the decoupling scale. Such a feature was also present in the Supersymmetry 
Parameter Analysis (SPA 05) prescription which was however based only on one-loop de- 
coupling and one-loop runninglfj The difference between the two results is larger than the 
experimental uncertainty by almost a factor of four. A peculiar feature of the presented 
results is that the MS-DR conversion and the DR running both depend on the unphysical 
coupling of the e-scalars to quarks. 

Higher order calculations are crucial in a variety of different contexts. However, at higher 
orders one has to deal with highly complex structures. Johannes Bluemlein reported on a 
method to simplify considerably physical quantities appearing in Quantum Field Theories 
(QFT) which depend on a single scale [5]. Single scale problems in massless or massive 
perturbative calculations in QFT can be expressed in terms of finite harmonic sums. These 
sums occur in the e expansion of the integrals for higher order corrections to QCD spit- 
ting functions and Wilson coefficients for space-like and time-like processes and high-energy 
scattering processes such as Bhabba scattering and many others. 

The idea is to perform a Mcllin transform and exploit the symmetry of Fcynman pa- 
rameter integrals in Mcllin space to perform simplifications. One derives algebraic and 
structural relations between multiple harmonic sums, and uses those relations to compactify 
considerably the result of higher order corrections. Applications of the method include the 
representation of 0(oq) Drell-Yan and Higgs-Boson production cross-sections, anomalous 
dimensions and Wilson coefficients at <D(ctg), asymptotic heavy flavor Wilson coefficients up 
to 0(al), and soft and virtual corrections to Bhabba scattering. Whether the method can 
be extended to multi-scale problems is something to explore in the future. 

At partonic level, a typical final state in small-x deep inelastic scattering off nuclei 
and hard proton-nucleus collisions can be characterized by the multiplicity of color-excited 
nucleons. Within the Reggeon field theory, each color-exited nucleoli is associated with 



a This is not consistent, one-loop matching needs a two-loop running. 
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the unitarity cut of the pomeron exchanged between projectile and nucleus. The pre-QCD 
Abramovsky, Gribov and Kancheli unitarity cutting rules formulated in the 1970s relate 
precisely these multi-pomeron exchange contributions to total, diffractive and inelastic cross 
sections. 

Koyla Nikolaev explained how, starting with an exact fc t -factorization for hard pQCD in 
nuclear environment, a dramatic revision of the AGK rules within QCD is found [6]. Two 
kinds of unitarity cut pomerons emerge which describe the color excitation and color rotation 
inelastic interactions. The departure of the unitarity cutting rules from the ones suggested 
by AGK stems from the coupled-channel features of non-Abelian intranuclear pQCD. In the 
Reggeon field theory language, the results entail a large variety of multi-pomeron couplings 
which vary from one universality class for hard pQCD processes to another. 

2.2 Jets, event shapes and unintegrated PDFs 

Jets enter a variety of measurements at colliders. What a jet roughly is is quite intuitive: 
a bunch of particles created in a hard interaction moving around a common direction and 
depositing energy on their way. Theoretically, jets are tools to create a link between the 
hadrons measured in detectors and perturbative calculations. Given a set of four-vectors, 
they are fully specified by the choice of a jet-clustering algorithm, its parameters and a 
recombination scheme. A variety of algorithms exist (the most simple classifications are 
cone types and sequential recombination type algorithms), some better than others as far as 
infrared safety, sensitivity to higher order or to hadronization effects are concerned. Once 
jets are defined, a number of observables can be studied. 

One interesting variable is the azimuthal correlation between two jets in QCD hard pro- 
cesses, A<p = \4>ji — <pj2\- In the presence of only two hard partons, without additional 
radiation, the two jets are back-to-back, A(f> = n, the non-zero azimuthal correlation pro- 
vides therefore a measurement of the additional QCD radiation. The azimuthal correlation 
is sensitive to a variety of physical effects: it depends on the jet algorithm and recombi- 
nation scheme; it is sensitive to soft/collincar gluon radiation and non-perturbativc effects; 
according to the recombination procedure in the clustering it can be affected by so-called 
non-global logarithms [7] and it is one of the observables commonly used for Monte Carlo 
(MC) tunings. Non-global logarithms have been the subject of many studies recently. The 
current status is that they can be resummed numerically in the large N c limit with single- 
log accuracy. Some time ago, it has been shown that the use of a jet-algorithm reduces 
the presence of undesired non-global logs [8] , however recently it has been shown that ad- 
ditional logs can originate from using a jet-clustering [9]. For the azimuthal correlation the 
recombination scheme turns out be to decisive to establish the presence of non-global logs: 
if one uses a recombination scheme which adds four-momenta, as is done at the Tevatron, 
the observable is non-global, if on the contrary one uses an E t weighted average azimuth, as 
was done at HERA, then the observable is continuously global. 

At last year's DIS workshop experimental results were presented by Magnus Hansson 
which did not agree with NLO 3-jet predictions [TU]. Triggered by this lack of agreement 
between data and theory theoretical work started to improve on the theoretical accuracy. 
The disagreement between fixed order perturbative calculations and data is to be expected in 
the back-to-back region since in this region multiple soft-collinear effects become important. 
Yazid Delenda reported on the resummation of soft-collinear logarithms in the azimuthal 



DIS 2007 



correlation distribution in the region n — A(j) ~ . Work in progress is now the matching 
with fixed NLO predictions, the inclusion of power corrections and the extension of this 
study to the Tevatron for which similar measurements exist. 



Event-shapes provide a continuous measure of energy and momentum of the final state. 
Usually, they are computed perturbatively at NLO and they have large soft-collinear loga- 
rithms which can be resummed at next-to-leading logarithmic (NLL) accuracy. Additionally, 
they have large hadronization corrections which give rise to 1/Q corrections, where Q is the 
hard scale of the process. These corrections can be modeled in terms of one universal pa- 
rameter ao- Universality implies that if one performs a combined fit of a s and ao then 
the values obtained should be independent of the event-shape used in the fit. This picture 
proved to work well both in e + e _ and DIS. However, it has been tested only in three-jet 
event shapes, those, like the thrust, whose first non zero contribution comes from events 
with three hard partons (one incoming and two outgoing in DIS). 

Andrea Banfi presented results for four-jet event shapes, whose first non-zero contribution 
is due to configurations with at least four partons in the event [12] . Specifically he considered 
-Kout, roughly the radiation out of the plane spanned by the three hardest partons and 
the D-parameter, proportional to the determinant of the energy-momentum tensor, 6 a ^ = 
J2hPhPh/(Q\Ph\)- Andrea presented fits of a s and ao from distributions of the /^-parameter 
for various values of y cu t, the cut used to select events containing three hard jets. It turns 
out that fits from the D-parameter give values of a s and ao which are compatible with other 
three-jet event-shapes, where the fit region is limited to the region to the right of the peak 
of the distribution. For K out on the other hand the shape of the distribution did not allow 
to perform a consistent fit. A possible explanation is that the K out distribution is sensitive 
to power corrections from the four-jet region currently not included. 

Final state predictions depend critically on which initial state partons enter the hard 
interaction and on the energy distribution of those partons. This is described by parton 
density functions. Standard PDFs depend on the factorization scale /i and on the longi- 
tudinal momentum fraction x, but they are integrated over the transverse k± component. 
Unintegrated parton densities (uPDFs) A{x, k\, (i) keep the information over the transverse 
component during the pcrturbative evolution. Cross sections are then given by the convolu- 
tion of uPDFs with off-shell partonic cross-sections a{xi,k\). If one considers e.g. the cross 
section for a single heavy flavoured parton as a function of pt, because of momentum con- 
servation at LO pt = 0, and non zero contributions appear first at NLO. If one instead uses 
uPDFs, the first non-zero contribution opens up already at LO. The advantage of uPDFs is 
therefore that by having the correct kinematics at LO one reduces considerably the size of 
the NLO corrections. At high energies the gluon density is dominant, therefore, as a first 
approximation, one can consider the contribution from gluons uPDFs only. 

The uPDF is determined by a convolution of a non-pcrturbativc starting distribution 
Aq(x, /io) and the CCFM evolution [13] A(x, k±,q) 




(1) 



where the distribution Ao is parametrized at a starting scale [Aq by 

xAo(x, no) = Nx~ Bs > (1 - x) 4 . 



(2) 
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Hanncs Jung presented results of fits of B g from inclusive Fi and from semi-inclusive i 7 ^ , giv- 
ing values olB g = 0.028 ±0.003 and B g = 0.286 ±0.002 respectively (for (j, = 1.2GeV) [1J]. 
The different value of B g from F£ changes the uPDF significantly, since the discrepancy is 
not covered by the experimental uncertainty. To resolve the discrepancy one can perform 
a fit to DIS dijet cross sections. The result for B g turns out to be very similar to the 
one obtained from F%. Regarding the residual discrepancy with the determination of the x 
component of A$ from F£ one can observe that the rise of the gluon distribution from F£ 
comes from lowest x points. If one considers only x > 2 ■ 10~ 4 and uses the F-i value of 
B g in the F^ distribution one obtains a x 2 /d.o.f. ~ 1.1. Finally, Hannes presented the first 
measurement of the intrinsic k t component of xA(x,/io) assuming a Gaussian distribution 
~ exp(— (k t0 — /j,) 2 /a 2 ). The fit to the intrinsic k t distribution, presented at this workshop 
for the first time, gives a ~ 1.5 and fi ~ 1.5 and is consistent with a Gaussian distribution, 
although other distributions are not excluded. One additional remark about uPDFs is that 
one can describe the A$ distribution between jets with Cascade using the new fit for the 
uPDFs. 

While NLO predictions for inclusive jet spectra typically have errors ~ 10 — 20%, NLO 
predictions for heavy jets^ fare much worse, the uncertainty being around ~ 40 — 60%. The 
NLO prediction for heavy-jets is thus a poorly constrained NLO calculation. Additionally, 
the experimental errors of the measurement at the Tcvatron arc smaller than the correspond- 
ing theoretical uncertainty. To understand the reason for the large theoretical uncertainty 
it is useful to examine the flavour production mechanisms. At LO the only production 
mechanism for heavy flavour is flavour creation (FC), 11 — > hh, where I is a generic light 
flavour and h denotes the heavy flavour under study. At NLO on the other hand two other 
processes enter, flavour excitation (FEX), Ih — > Ih, where the incoming heavy quark stems 
from a collincar splitting of an incoming gluon, and gluon splitting (GSP) 11 — > 1(1 — ► hh). 
Since the NLO processes are enhanced by collinear and soft logarithms, ln(p t .j C t/™Q) (where 
rriQ denotes the mass of the heavy quark), they turn out to be more important that the LO 
contribution. This means that a NLO calculation, which treats FC at NLO, but where FEX 
and GSP enter only at LO has very large K-factors (~ 5 at the Tevatron and ~ 8 at the 
LHC, for accessible p t values) and therefore also very large scale dependencies. 

Giulia Zanderighi showed that if one instead clusters events using an infrared-safe, flavour 
algorithm, which takes into account the difference between quark and gluon QCD pro- 
duction, then flavour becomes an infrared safe quantity. This implies that no logarithms 
ln(pt.jet/wQ) are present and one can accordingly perform a light-quark calculation e.g. with 
NLOJET+- 1- [15] . This gains roughly a factor 3 in accuracy: K factors become ~ 1.3 and 
scale variation uncertainties are around 10 — 20%, as is to be expected from a true NLO 
calculation [16j . The method however requires a good understanding of single and double 
6-tagging efficiencies, and further experimental investigation in this direction is needed. 

2.3 Parton showers and matrix element event generators 

QCD Monte Carlo (MC) generators and simulations are vital for physics at the LHC. It 
is therefore crucial to examine critically different MC components, including perturbative 
aspects (parton showers, matrix element corrections, matching), non-perturbative aspects 

b Experimentaly a heavy-jet is any jet containing at least one heavy tagged object. 
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(hadronization, underlying event) and tuning of the event generators. Parton showers (PS) 
reflect our understanding of pQCD to all-orders. They are commonly believed to capture 
at least the leading logarithmic structure. Observables sensitive to radiation in a limited 
phase space region (energy flow distributions, event shapes, . . . ) are sensitive to non-global 
logarithms even at leading order. Is it therefore important to re-examine the accuracy of 
showers in these instances in order to establish if for example leading non-global logarithms 
are erroneously tuned in the MC parameters. 

The starting point of the study presented by Mrinal Dasgupta is the fact that angular or- 
dering (AO) catches the relevant part of these non-global logarithmic corrections. HERWIG, 
therefore, which is based on angular ordering is expected to be close to the full resummation 
of the non-global logarithms (known only in the large N c limit). PYTHIA (before v. 6. 4) 
uses as an ordering variable, the virtuality m 2 and rejects a posteriori configurations which 
do not respect AO, and is expected to do worse. PYTHIA v. 6. 4 on the other hand, like 
ARIADNE, is based on dipole phase space, and should have the full leading logarithmic 
behaviour. In the study presented by Mrinal Dasgupta these expectations are tested by 
comparing full resummed results matched to NLO with predictions for parton showers for 
the radiation in a slice in rapidity Ar) as a function of E t , the total transverse energy of 
hadrons in the rapidity slice [17] . It turns out that the above expectations are confirmed 
unless one takes a quite large region of pseudorapidity, e.g. Arj = 3. For large rapidity slices 
a large discrepancy is found between the new PYTHIA (v. 6. 4) and the resummed result. 
Further studies are needed to clarify this discrepancy and to generally understand the vari- 
ous parton showers in a quantitative fashion. In general, it has been stressed that whenever 
possible one should compare results from HERWIG and PYTHIA, and when possible of 
other Monte Carlos. 

As stressed above, Monte Carlos today are indispensable tools for many experimental 
studies. One of their very conventional application is to estimate non-perturbative (NP) 
effects. NP effects at hadron colliders include both final state hadronization effects and 
effects due to the underlying event (UE). Single jet inclusive distributions, parametrized 
in terms of a jet-radius R = \J (Ay) 2 + (A</>) 2 have l/ptjct power corrections due to the 
hadronization and UE. 

Lorenzo Magnea [TS] explained that it is possible to show analytically that hadronization 
corrections are distinguishable from UE because they exhibit a singular dependence on R, 
while UE effects are proportional to the radiation they collect in a given region and so scale 
as ~ irR 2 . Lorenzo then presented a MC study where one runs a MC at parton level (p), 
hadron level without UE (h) and with UE (u). One then selects events with two hard jets 
with the hardest jet in a chosen p t range and one defines for each hadron level a measure of 
the size of the power correction 

A^ /U) = \ ( P ( T h { u) + V% U) - &l ~ P { £) , AjT h) = Af£ J - (3) 

One can examine Ap^ divided by the leading behaviour in R (1/R for hadronization 
effects and R 2 for UE) for two different partonic channels, e.g. a quark dominated channel, 
qq — > qq, and a gluon dominated one, gg — > gg, both for HERWIG and PYTHIA. With 
HERWIG the result of this operation is a flat distribution as a function of R, signaling 
that the leading R dependence has been correctly extracted. Additionally, one can see that 
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the size of the UE is independent on the hard partonic scattering channel. With PYTHIA 
on the contrary the UE seems to "know about the hard scattering" i.e. instead of a flat 
distribution one has a leftover l/R dependence. This is perhaps surprising. Additionally the 
UE is sensitive to the hard channel (increases when going from quarks to gluons). A correct 
modeling of the UE is the basis for precision phenomenology at the LHC, since for example 
this input enters into the determination of the jet energy scale. These types of discrepancies 
should therefore be addressed soon. The work presented here is based on disentangling 
hadronization from UE effects by exploiting their different R dependences. It was therefore 
recommended that, whenever possible, measurements for different values of the jet-radius R 
should be provided. 

In a more general context parton showers rely on the universal soft and collinear (SC) 
factorization of the QCD matrix elements. This is a universal property and is true to 
all perturbative orders. In an ideal world this should be the only approximation made 
in MCs. However, practically all current MCs are subject to many other approximations, 
e.g. interference diagrams are treated only approximately with angular ordering, the color 
algebra is valid only in the large N c limit, spin treatment is usually inexact, other first 
approximations arc carried in order to simplify phase space and further arbitrary techniques 
are often employed. 

Zoltan Nagy reported on ongoing work to formulate a parton shower which removes 
all the above approximations [19] . The method is based on recursive equations that can be 
used to generate a lowest order parton shower for hard scattering in hadron-hadron collisions. 
The formalism is based on the factorization of soft and collinear interactions from harder 
interactions. It incorporates quantum interference between different amplitudes in cases 
where the interference diagrams have leading soft or collinear singularities. It incorporates 
the color and spin information for the hard partons. One of the motivations for this more 
formal development of a parton shower is to have a method that can be merged to NLO 
calculations in a natural way. 

While parton showers are crucial for many studies, it is also well-known that in many 
cases parton showers fail dramatically. This is for instance the case in events with many hard 
jets: a parton shower has the full matrix element only for the primary hard scattering involv- 
ing few partons, while all other emissions are treated in the soft-collinear approximation^ 
Therefore, any study which looks at properties of events with many hard well-separated 
particles in the final state should be based on tools which go beyond the soft /collinear 
approximation of parton showers. Today tree level matrix element generators (ALPGEN, 
Helac, MadGraph/MadEvent, Amegic++, . . . ) exist which can treat according to the pro- 
cess up to around 5-8 particles in the final state including full spin correlations/interference. 
Simon Visscher reported on progress in Madgraph/MadEvent (MG/ME) [20]. Specifically, 
MG/ME have now two different matching schemes, an MLM type based on event rejection 
after the PS and a CKKW type based on Sudakov rewcighting and vetoing on showers. 
This work set a basis for a systematic comparison between different generators, matching 
techniques and shower algorithms. |21j 



c Traditionally, using the soft /collinear approximation means that there is a deficit of radiation in the 
large angle region. This is however not always the case, one can enhance the shower arbitrarily so as to fill, 
even overfill, those regions, but this is artificial and ad-hoc. 
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As an example Simon presented results obtained with MG/ME for ti + 0, 1, 2,3 . . . jet 
events. These processes are important for a variety of reasons. In addition to their interest 
per se, they are a background to ttH production and in two-Higgs Doublet Models (2HDM) 
pp — > W + W~bbbb could be the most interesting channel to discover the charged Higgs. In 
this case one needs a reliable ti + 0, 1, 2, 3 . . . jet event sample. If one looks at the rapidity 
distribution of the hardest jet, predictions obtained with Pythia are depleted of radiation 
in the central rapidity region. Full matrix element based predictions are not, showing the 
need to go beyond SC approximation of parton showers for this type of analysis. 

Despite the fact that QCD and EW corrections are nothing but higher order corrections 
dictated by the same SM Lagrangian, traditionally the two corrections have mostly been 
treated separately Alcssandro Vicini reported on a calculation of combined QCD (0(a s ) + 
PS)) and EW (C(a ow ) + PS) corrections to charged current Drell-Yan (DY) processes at 
hadron colliders [52]. DY processes are unique at hadron colliders. The presence of high p t 
leptons, together with the large cross section makes their detection very easy. DY processes 
are useful for validating pdfs, as luminosity monitors or to obtain precision measurement 
of EW observables (Mw ■■■)■ W/Z production in association with jets is an important 
background to SM and beyond Standard Model (BSM) signals, especially to new gauge 
bosons. In pure QCD the state-of-the-art are fully differential NNLO predictions exist [23] . 

Alessandro Vicini presented results obtained with the HORACE event generator, which 
includes state of the art EW radiative corrections to DY processes including exact 0(a ew ) 
radiative corrections matched to multiple photon radiation via PS. It's a fully exclusive event 
generator, which can be easily interfaced to QCD showering programs like HERWIG, events 
are saved in a Les Houches compliant form and it can be interfaced to the LHAPDF package. 
The QCD and EW corrections are matched with an approximate additive procedure, which 
works to 0(a s a ew ) as long as hard non-collinear radiation effects are not important. Beyond 
the additive approximation a full two-loop 0(a s a cw ) calculation is needed. Results presented 
for Mjf and p l ± show that around the jacobian peak positive QCD corrections compensate 
negative EW corrections and that EW effects arc mandatory in order to extract Mw, only 
QCD PS are not sufficient. Another interesting QCD-EW interplay effect lies in the fact the 
convolution with QCD PS modifies the relative size and shape of the EW corrections. Due 
to the presence of EW logarithms, EW effects become important in the tails of distributions 
(My ~ 2TeV, p^_ ~ ITeV) they can amount to up to ~ O(30) % corrections, however 
cross-sections are very tiny in those regions. 

3 Particle production in ee, ep, pp and pp collisions 
3.1 Pentaquark searches 

In 2004 a comprehensive program to search for pentaquarks (PQ) in high statistic, high 
resolution experiments was started by the CLAS collaboration at Jefferson Lab. P. Rossi 
(CLAS) presented recent results for the four photoproduction channels studied so far [24] . 
No evidence for a PQ signal was observed. From these results CLAS set an upper limit on 
the PQ production cross section on protons and neutrons. 

M. Del Degan (HI) presented the invariant mass spectrum of the S-7T system, studied 
using deep inelastic scattering (DIS) data collected with the HI detector at HERA [25]. A 
clear signal for the well established unexotic S(1530)° baryon is observed. Despite having 
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similar statistics as NA49, the S~~/° PQ signal could not be confirmed. Upper limits 
for PQ production from HI are in agreement with those published by the ZEUS 

collaboration. 

3.2 Deuteron and antideuteron production 

A first measurement of the production of deuterons and antideuterons in DIS was presented 
by S. Chekanov (ZEUS) [26 . It is interesting to mention that this is the first measurement of 
deuteron production in elementary particles collisions. Production rates of (anti)deuterons 
are significantly reduced relative to the production of (anti)protons, consistent with other 
world measurements. The production rate of deuterons was measured to be approximately 
three times larger than for antideuterons. The production rates were also studied in terms 
of the coalescence model. 

3.3 Exclusive final states in e + e 

BaBar has now taken over 400/6~ 1 of data and has used this data in a variety of tests 
of QCD using different production channels. The first observation of e + e~ — > p p and 
e + e~ — > cf>p was reported by S. Saremi (BaBar) [27] , The final states for these channels 
are even under charge conjugation. Such processes can be understood in terms the Two- 
Virtual-Photon Annihilation model. Other results presented were a new test of factorization 
in the BO — » D* + ojit~ channel, for which there is good agreement, and a measurement of 
the r\ and r( transition form factors, were the ratio of the form factors is inconsistent with 
theoretical predictions. 

3.4 Identified particle measurements in ep and pp 

Recent results on strange particle production (K ± , Kg, A's) and Bose-Einstein correlation's 
(BEC) between kaons were presented in a talk given by B. Levchenko (ZEUS) 2SJ . The 
parameters for the BEC agree well for different kaons and are consistent with HI and e + e~ 
results. No sizeable barion-antibarion asymmetry was observed. The ratio of baryons to 
mesons in the resolved photoproduction regime is much larger than in e + e~ and is not 
described by Pythia predictions. 

This observation agrees well with similar results from pp, presented by M. Heinz (STAR) . 
STAR also presented the px spectra for different particles, which for the first time can be 
described by the NLO pQCD predictions due to improved fragmentation functions (FF) 
for baryons and strange particles using the light-flavour separated measurements in e + e~ 
collisions from OPAL. 

3.5 Fragmentation functions at HERA 

The scaled momentum spectra in the current region of the Breit frame in DIS ep scattering 
at HERA was measured by ZEUS, presented by B. Brzozowska, with high precision data 
covering a large range of energy scale (5 to 170 GeV) [30]. The NLO pQCD theoretical 
predictions cannot reproduce the data in the entire phase space and cannot describe the Q 2 
evolution of the xp distribution. Differences between predictions using different FF were 
found to be small. 
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3.6 Particle production in jets 

The relative roles of perturbative and non-perturbative QCD in the development of jets 
have been studied by CDF. In his talk L. Pinera presented the momentum distributions and 
multiplicity of charged particles in jets, momentum correlations of particles and their kx 
distribution [31]. The energy scale range covered in this analysis goes from 20 to 160 GeV. 
The data are generally well described by pQCD within MLLA approach, but at large kr 
the theory predictions significantly deviates from data. 

4 Particle production in collisions with heavy-ions 

4.1 Property of the Quark-Gluon Plasma 

Hadronic final states and their correlations were used by the PHENIX collaboration (talk 
presented by C. Ogilvic) to investigate properties of the Quark-Gluon Plasma (QGP) and 
it's impact on particles passing through it |32j . Energy loss and meson suppression deliver 
information on the density of QPG while direct photons are insensitive to the QPG. These 
properties can aid the comparison of collisions of protons with that of heavy-ions. The 
broadening of jets can be explained in terms of induced gluon radiation in the QGP. High 
Pt jet production is biased towards being produced on the surface of the QGP, the second 
jet in these types of events has to then travel through the QGP to the far side. Studies 
of the far side reveal that the second jet is suppressed and takes on the properties of the 
medium it travels through. Studies show that the medium also responds to the passage of 
the second jet and heavy quarks are seen to lose less energy than light quarks. 

4.2 Colour transparency 

A study of colour transparency (CT) was presented by M. Holtrop (CLAS) [33]. The CT was 
studied in terms of nuclear transparency T(A,Q 2 ) and is clearly observed in p° electropro- 
duction, where T(A, Q 2 ) was measured to increase with increasing Q 2 , in a good agreement 
with predictions of the theoretical model by Kopclovich et al. 

4.3 Hadronisation in nuclear environment 

Recent results from HERMES were presented by Z. Akopov and Y. van Haarlem [34l [35] . 
Detailed studies of hadronisation on many nuclear targets were performed. Substantial 
nuclear attenuation was observed as a function of different kinematic variables. For the 
first time these studies were done using double-differential distributions. For the first time 
the formation length of the nuclear attenuation was studied. Absolute measurement of the 
PT-broadening has been performed using different targets and hadron types for different 
kinematic variables. 

5 Jet production 
5.1 Prompt photons 

S. Chekanov (ZEUS) presented recent results on prompt photons with an associated jet 
in photoproduction at HERA [36] . Prompt photon identification was based on conversion 
probabilities measured by a dedicated detector - Barrel Pre-samplcr. Both PYTHIA and 
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HERWIG fail to describe the shape and normalization of the cross sections. Best agreement 
is achieved using the /c^-factorization approach. An improved agreement with other NLO 
pQCD calculations was achieved by hardening the cut on the transverse momentum of the 
photon. 

The HI collaboration (talk by K. Miiller) has performed prompt photon studies in DIS 
[37j . The Prompt photon signal was extracted using a log- likelihood approach based on a 
shower shape analysis. The measurements were done both for isolated photons and photons 
with jets. The LO MC models underestimate the data but the shapes are generally described. 
NLO pQCD calculations were not yet available. 

O. Atramentov (D0 ) presented inclusive photon+jct results in pp collisions [35]. Triple 
differential cross sections were measured as well as the ratios of these cross sections, sig- 
nificantly reducing the effect of experimental and theoretical correlated uncertainties. Four 
different photon+jet topologies were studied to explore different kinematic regions of the 
gluon distribution functions. For some kinematic regions quantitative deviations from cur- 
rent QCD predictions are observed. 

5.2 Jets photoproduction at HERA 

High-i^T jets in photoproduction (PHP) measured by ZEUS and presented by H. Perrey show 
sensitivity to the gluon PDF and have the potential to further constrain the parton densities 
of the proton and photon [39] . None of the compared photon PDFs provided an adequate 
description of the ZEUS resolved data, although the direct enriched cross sections, which 
are less sensitive to the photon PDF, are in good agreement with NLO pQCD calculations. 

The photoproduction of dijet events with a large rapidity gap between the jets was used 
by ZEUS to estimate the colour-singlet exchange contribution, as reported by A. Savin 
(ZEUS) [40]. Multi-jet photoproduction was used to study the properties of Multi-Parton 
Interactions (MPI) at HERA. The low-invariant-jet-mass region can only be described by 
adding the MPI to the MC simulation, thus providing a good testing ground for different MPI 
models. This information is very important for understanding the structure of underlying 
events at the LHC. 

5.3 Study of underlying events in DIS 

Further studies on underlying events at HERA were performed by HI using the DIS data 
by measuring jets with low transverse momenta. S. Osman (HI) demonstrated in his talk, 
that the standard QCD MC models fail to describe the data in different azimuthal regions 
with respect to the leading jet [41] . Adding the MPI significantly improves the description 
at low Q 2 where resolved photon processes are expected to contribute, but still fails at high 
Q 2 . Due to lack of NLO pQCD calculation suitable for this analysis no conclusion from 
high-order contribution was drawn. 

5.4 Charge current jets at HERA 

The first charge current (CC) ep jet analysis of the HERAII data was presented by H. Wolfe 
(ZEUS) [42]. For the first time the three- and four-jet events in CC DIS are analyzed. 
Inclusive jets and multi-jet cross sections were measured. Total polarized inclusive jet cross 
sections agree with predictions of the Standard Model. 
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5.5 Jets in ee 



An investigation of high momentum hadron and jet production by the OPAL collaboration 
was presented by A. Krasznahorkay [43] . Previously measured exclusive jet cross sections 
by L3 was not well reproduced by the pQCD NLO predictions. The new measurement by 
OPAL is well described by the theory. When compared in the same kinematic regime as L3, 
the OPAL and L3 measurement do not agree well and this difference has to be understood. 
The pt spectrum of hadrons is well described by the calculation except in the very high pr 
bin. 

5.6 Jets in pp and pp 

A preliminary measurement of the bb dijct production cross section and the angular correla- 
tion of the jets was presented by S. Vallecorsa (CDF) [33] . The NLO calculations describe the 
data well. The bb angular correlation shows that events are mainly produced by the flavour 
creation mechanism. The low A(f> tail suggests non-negligible contribution from other pro- 
cesses. It was demonstrated that the inclusion of the underlying event in the simulation 
significantly improves the description of the data. 

M. D'Onofrio (CDF) reported on measurements of VF+jets and Z+jets production cross 
sections at CDF [35]. Boson+jet production channels are fundamental to test pQCD, un- 
derlying events and new theoretical calculations. Comparisons to different calculations were 
shown. The size of non-perturbative corrections was also estimated. The measurement is 
important since it gives an estimate of background for new physics like top, Higgs and SUSY 
searches. 

Jet production measurements was also reported by J. Cammin (D0 ) [46] . Good agree- 
ment with NLO pQCD predictions over a large px range was demonstrated. Systematic 
uncertainties are dominated by the jet energy scale. Measurements of jets with a \i tag, 
used to identify Heavy Flavor jets, and Z+jet production were also presented. 

Prospects for inclusive jet cross section measurements with early data at ATLAS were 
presented by D. Clements [47]. Experimental errors are expected to be dominated by the jet 
energy scale. Theoretical errors at high px are dominated by the uncertainty on the high-x 
gluon PDF. 

5.7 Inclusive jet production and jet correlations in DIS 

Recent results on angular correlations in three-jet events and jet substructure were presented 
by E. Ron (ZEUS) [35]. Fixed-order calculations, separated according to the colour config- 
urations, were used to study the sensitivity of the angular distributions to the underlying 
gauge structure. Using different variables one can distinguish between contributions from 
triple-gluon vertexes in quark-induced and gluon-induced processes. The measurements are 
consistent with the admixture of colour configurations as predicted by SU(3) and disfavour 
some other combinations. Subjet structure is reasonably described by the fixed-order QCD 
calculations and are consistent with the dominance of quark-induced processes. 

Inclusive jet production in DIS at high Q 2 is a well established measurement and is well 
reproduced by different pQCD calculations. This is why this regime is used for precise QCD 
measurements and extraction of as- An updated HI analysis with improved precision was 
presented by T. Kluge (HI) [49] , The as fit was performed minimising the experimental 
uncertainty. A first extraction of as using the ratio of inclusive jet to the inclusive DIS 
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cross sections was made which provides a more precise measurement. The ZEUS measure- 
ments were presented by T. Schoerner-Sadenius [50] . Together with inclusive cross sections 
ZEUS measured dependence of inclusive jet cross section on jet radius and used this mea- 
surement to extract the a s value. The inclusive dijet sample was analyzed using combined 
HERAI and HERAII data, correspondent to integrated luminosity of 210 pb~l. For all 
these measurements the theoretical uncertainties still dominate. 

Inclusive jet production cross section in the low Q 2 region, 5 < Q 2 < 100 GeV 2 , was 
presented by A.Baghdasaryan (HI) [5T]. The analysis extends to lower Q the similar analysis 
performed by HI at high Q 2 . Good agreement with NLO pQCD calculation was observed for 
Q 2 above 10 GeV 2 . For further talks on jet production in the low- a; regime see forward jet 
production and HERA multi-jet production at low- a; and low-Q 2 presented during a special 
Structure Function WG session. 

6 Summary 

In addition to several formal theoretical developments, further experimental and theoretical 
contributions demonstrated great progress towards a more precise understanding of QCD 
final state production. Enormous progress in the theoretical calculations triggered by more 
precise data included a first resummation for QCD 3-jet production at Next-to-Leading 
Order which should provide a better description of dijet azimuthal correlation measurements 
both at HERA and the Tevatron. An additional dedicated session was held to discuss 
as, where theoretical uncertainties dominate experimental measurements at HERA. NNLO 
calculations for 3-jet production in e + e~ were presented and are expected to lead to new 
fits for a s with reduced theoretical uncertainties. 

Future new physics searches at the LHC will rely on the most precise determinations of 
the parton density functions. HERA and Tevatron photon + jet production measurements, 
highlight regions where continuing improvements to theoretical models and inputs are still 
required and measurements of high Et forward dijet production at HERA and inclusive 
jet production at the Tevatron should also help further constrain these inputs. Progress 
towards a precision determination of unintegrated PDFs and a first determination of their 
intrinsic kx component was also presented at the meeting. 

Tests of power corrections arising from hadronization effects have been extended to 3-jet 
event shape distributions. Progress in perturbative evolution and hadronization issues, the 
role of the underlying event implementation in Monte Carlos, reducing model approxima- 
tions for parton showers and issues that arise when Monte Carlos are tuned to non-global 
parameters were highlighted. The importance of multiple interaction modeling was demon- 
strated by HERA 7p multijet production measurements and bb jet azimuthal correlations 
at the Tevatron. Improvements using a new flavor jet algorithm to theoretical uncertainties 
for heavy flavor jet production were demonstrated. 

Recent exclusive final state measurements include an observation of two virtual photon 
annihilation final states at BaBar. Work in progress to resolve pentaquark issues currently 
show no further evidence for pentaquark states. A first measurement of the ratio of deuteron 
and antideuteron production rates in DIS at HERA is smaller than expected in the coales- 
cence model. NLO pQCD predictions using improved fragmentation functions are found 
to describe the measured particle p t spectra at RHIC, however the effect of fragmentation 
functions is small and unable to account for observed differences to NLO pQCD predictions 
in scaled momentum spectra measurements in DIS at HERA. In collisions with heavy ions 
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recent measurements at RHIC on the impact of Quark-Gluon plasma properties to particles 
passing through it were presented, together with first double differential nuclear attenuation 
measurements on a series of different nuclear targets at HERMES and color transparency 
effects which were reported in p° electroproduction measurements at CLAS. Further details 
for each individually presented topic can be found in the dedicated contributions to these 
proceedings. 
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